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Abstract 

The particle-antiparticle asymmetries of A production in 250 GeV/c vr^, K^, and p -nucleon 
cohisions are studied with two model parametrizations of quark to A fragmentation functions. It is 
shown that the available data can be qualitatively explained by the calculated results in both the 
quark-diquark model and a pQCD based analysis of fragmentation functions. The differences in 
the two model predictions are significant for beams, and high precision measurements of the 
asymmetries with detailed xp and Pt information can discriminate between different predictions. 
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Much of the recent interest in the structure of the A hyperon stems from the fact that 
its quark structure is rather simple in the naive constituent quark model, since it contains 
three different quark flavors {u, d, s), and furthermore its spin is carried only by the strange 
quark. Therefore any departure from this picture signals the presence of relativistic and non- 
perturbative effects. Experimentally the most effective means to investigate the A structure 
is through the quark to A fragmentation in various processes 

Unfortunately there are still no available measurements about the relations between differ- 
ent flavor-dependent quark to A fragmentation functions, such as the relation between the 
favored and unfavored fragmentation functions, and the relation between the flavor struc- 
ture of the favored fragmentation functions D^{z) and D^[z). Nevertheless there has been 
recent new measurement of the A-A production asymmetries in hadron-nucleon collisions 
with several different hadron beams at 250 GeV/c by the Fermilab E769 Collaboration j^. 
The purpose of this work is to show that the E769 data of A-A asymmetries in hadron- 
nucleon collisions with different hadron beams may provide useful information about the 
flavor structure of quark to A fragmentation functions. 

We first look at the formalism for inclusive production of a hadron C from hadron A and 
hadron B collision process 

A + B^C + X, (1) 
where the Mandelstam variables s, t, and u are written as 

s = {Pa + Pb? = Ml + M| + 2Pa ■ Pb, (2) 
t = {PA- Pc? = Ml + Ml- 2Pa ■ Pc, (3) 
u={Pb- Pc? = Ml + Ml- 2Pb ■ Pc- (4) 

The experimental cross sections are usually expressed in terms of the experimental variables 
xp = 2Pi^j ^ and Pt at a given s, where Pl = xf\/s/2 and Pt are the longitudinal 



and transversal momentum of the produced hadron C with energy Ec = \/ -|- P^ = 
a/ + Pi + P^ = a/ M^, + Pr^ + x|.s/4 in the center of mass frame of the collision process, 
and s is the squared center of mass energy. In the coUinear factorization theorem, the 
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kinematics of the subprocess a + b^c + d is defined as 



Pa = XaPA, (5) 

Pb = XbPs, (6) 
Pc = Pc/z, (7) 

from which we get the parton level Mandelstam variables s, i, and u 

S= {Pa+ Pbf = '^Pa ■ Pb = XaXbS, (8) 
i = {Pa - Pc? = -2pa ■ Pc = ft, (9) 

u = {Pb- Pc? = -2pb ■ Pc = fu. (10) 
The cross section can be written in terms of the parton level subprocess as 

da=^—- = / 2 fai^a)fb''iXb)s6{s + t + u)—^{Xa,Xb,z)D^{z), 

a,b,c,a 

(11) 

where fai^a) and f^{xb) are the parton distributions in the beam hadron A and target 
proton B with Bjorken variables Xa and Xb respectively, D^{z) is the fragmentation function 
of parton c into the produced hadron C with energy fraction z of hadron C relative to the 
scattered quark c, and ^^-^f—{xa, Xb, z) is the subprocess cross section. From the 5{s + i + u) 
function, we get the energy fraction z 

Xgt -\- XbU (12) 



XaXb^ 

and the two Bjorken variables for quarks a and b in hadrons A and B 

XbU Xcit /io\ 
Xa = , Xb = . (13) 

XbZS + t XaZS + U 

The integration over Xa and Xb can be done for both of Xa and Xb in the range [0, 1], under 
the constraint condition < z < 1, or equivalently, for Xa and Xb in the ranges [x™", 1] and 
[x™", 1] respectively, with 

a:r" = -^^, a;r" = -^. (14) 
Xb-s + r s + u 

The hard scattering framework of Eq. (fTT|) is perturbative QCD based and therefore 
strictly valid only for high Pt > 1 GeV. Nevertheless, it has been recently argued that QCD 
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factorization is also applicable for semi-inclusive processes at low transverse momentum , 
therefore we make extrapolations to smaller Pt in this paper for the purpose of illustration. 
The kinematical expressions that we use, however, are different from the conventional ex- 
pressions in terms of rapidity ll| and are valid even for small Pt- 

The cross sections of the subprocess at the parton level, a{s,i,u), can be found 
in Refs. We adopt the Gliick, Reya, and Vogt (GRV) leading order unpolarized 

parametrization for the nucleon parton distributions ^13|. For the parton distributions of 
the and vr", we employ the parametrization forms of Refs. and ^|, respectively. 
The above quantities are well constrained by a vast number of available experimental data, 
so we can focus our attention on the less known quark to A fragmentation functions. 

We parametrize the quark to A fragmentation functions D^{z) by adopting the Gribov- 
Lipatov relation Q| 

D^iz)^q\x), (15) 

in order to connect the fragmentation functions with the quark distribution functions q^{x) 
of the A. More explicitly, we adopt a general form to relate fragmentation and distribution 
functions, as follows [17[ 

Dliz) = Csiz)z-qliz), 
where a distinction between the valence (V^) and the sea (S) quarks is explicit. The above 
formulae are always correct, since Cv{z) and Cs{z) are in principle arbitrary functions. 
We should consider Eq. as a phenomenological parametrization for the fragmentation 
functions of quarks and antiquarks, as follows 

D^iz) = D^iz) + Dliz), 

Three options were found ^y}\ to fit quite well the available experimental data of proton 
production in e~^e~ inelastic annihilation: (1) Cy = 1 and Cs = for a = 0; (2) Cy = Cs = 1 
for a = 0.5; and (3) Cy = 1 and C5 = 3 for a = 1. We adopt these three options to reflect 
the relation between unfavored and favored fragmentation functions of the A. 

There is no direct measurement of the quark distributions of the A. But we can relate the 
quark distributions between the proton and the A by assuming SU(3) symmetry between 
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the proton and the A 



for valence quarks, and 



(18) 

y{x) = |mv'(x) - 



u^{x) = (/^{x) = I [u{x) + s{x)] 



(19) 

s^(x) = d(x), 

n 

for sea quarks. We adopt the GRV parametrization [13| of the quark distributions q{x) of 
the nucleon. In this way, we get a complete set of quark distributions in the A with both 
valence and sea quark distributions. 

It is well known that the flavor structure of u and d quark distributions of the proton is 
diffe... between the ,ua..U-di,ua.. .ode. UUU and a pQCD based analysis QQ: 
the quark-diquark model predicts that d{x)/u{x) — > at x ^ 1 whereas a pQCD based 
approach predicts that d{x)/u{x) ^1/5. A discrimination between the two models requires 
very high precision measurement of the structure functions at large x and is difficult. On the 
other hand, it has been also shown [2J that this flavor structure of the quark distributions 
at large x is even more signiflcant in the case of the A, with a large difference between the 
ratio of u^{x)/ s^{x): the quark-diquark model predicts that u^{x)/ s^{x) — > at x — 1, 
whereas the pQCD based approach predicts that u^{x)/ s^{x) 1/2. This will produce a 
large difference in the ratio of fragmentation functions D^{z) / D^{z)^ which might be more 
easily accessible experimentally via quark to A fragmentation 24| . 

The valence quark distributions of the A in the quark-diquark model and the pQCD 
based analysis have been explicitly studied 

ESQ 

and we adopt the parametrizations 
given in Ref. [2^. To describe the A fragmentation, it is important to take into account 
the sea contributions in the model construction. In order to use the sea quark distributions 
from other parametrization while still keeping the flavor structure of the valence quarks 
as predicted in the two models, we re-scale the valence quark distributions by a factor of 
su{3)(-^)/'^v,thi^)^ where the subscript "^^/(S)" denotes the valence quark distributions 
of the A in the SU(3) symmetry model |3| and "t/i" denotes the corresponding quantities 
predicted in the quark-diquark model or the pQCD based analysis {2^. This is done in order 
to normalize the A quark distributions to well known proton quark distribution parametriza- 
tions. Notice that the valence w-quark distribution then becomes that of the SU(3) model, 
while the others get a rescaling factor. In this way we can adopt the sea quark distributions 



from the SU(3) symmetry model as the sea distributions in the quark-diquark model and the 
pQCD based analysis, to reflect the contribution from the unfavored fragmentation. Thus 
we get another two sets of quark to A fragmentation functions, denoted as the quark-diquark 
model and the pQCD based analysis later on. This procedure is done with the main moti- 
vation of constructing realistic quark to A phenomenological fragmentation functions, which 
has some features coming from specific theoretical arguments, i.e., the quark-diquark model 
and the pQCD based analysis. 

With all of the above mentioned subprocess cross sections at the parton level, parton 
distributions for both the beam and the target, and also the quark to A fragmentation 
functions, we can calculate the x^-dependent A- A asymmetries in hadron-nucleon collisions 

pmax 

A{xf) = / " d'Pr 

J pmin 

and those with P^-dependence 
A{P^) = I dxp 

We do not introduce any additional parameters to fit the data for the input quantities. 

The E769 Collaboration measured the asymmetries with 250 GeV/c vr^, K^, and p 
beams on the proton target. The data are expressed as functions oi xp and over the 
ranges —0.12 < xp < 0.12 and < P^ < 3 (GeV/c)^ for positively charged beam, and 
—0.12 <xp< 0.4 and < P^ < 10 (GeV/c)^ for negatively charged beam. In fact, the xp- 
dependence was measured without Pp cut, and the P^-dependence was measured without 
Xp cut. Notice that in the Pp integration we have considered the small (< 1 GeV) region, 
where the hard scattering formalism is not strictly valid. We have checked that the final 
results are not really sensitive to this lower limit. 

In the theoretical calculation, we need to set the limits of P^*" and P^"^ for the xp- 
dependence, and of x^*" and x^"^ for the P^-dependence. We will choose the above men- 
tioned ranges in the calculation. In Fig. HI we plot the calculated results for the XF-dependent 
asymmetries A{xp\ and find that both the quark-diquark model and the pQCD based anal- 
ysis can explain the trend of the data. The quark-diquark model seems to be somewhat 
more favored by the data for beams, especially at large xp. The differences between 
the quark-diquark model and the pQCD based analysis can be understood by the fact that 



dV^ 



d^ 



d^p 



d3 



3^A 



T 



+ 



dV 



pmin 



d'^PA d'^p 



(20) 



d3 



cr 



d3 



d^pA d^p^ 



dxi 



dV^^ dV 

d^pA d^p^ 



(21) 
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xf ^ Xa — Xb in the Bjorken limit, and that the large xp behaviors is mainly controlled 
by the s and s quark fragmentation of the beam at large xp- In the quark-diquark model, 
the s quark fragmentation dominates over the light-flavor quark fragmentation. Therefore 
the dominant s fragmentation of the beam favors A production, whereas the dominant 
s fragmentation of the beam favors A production. This can explain why the asym- 
metries A{xf) for the K~ beam is negative at large xp whereas they are positive for the 

beam. From the different flavor structure of D^[z) / D^[z) ^ in the quark-diquark 
model and of D^{z)/ D^{z) 1/2 in the pQCD based analysis, we can explain why the 
quark-diquark model predicts large magnitude for the asymmetries A{xf) at large xp for 

beams, since the s quark fragmentation dominates over the light-flavor fragmentation 
in the quark-diquark model. We thus conclude that the large xp behavior of the A- A asym- 
metries can discriminate between the two model predictions. We present in Fig. |21 of the 
calculated results of the A{xp) asymmetries with two different options of P^-cut: the inte- 
grated ranges of [(P™"^)^, {P^^'^Y] are [0,5] and [5, 10] respectively. We find that the trend 
in the XF-dependence does not change much for the two options, although the magnitudes 
are different. 

We also plot the calculated results for the P^-dependent asymmetries A{P^) in Fig. 13 
and find that both the quark-diquark model and the pQCD based analysis can explain the 
P^-dependence trend, but the results differ slightly in magnitude from the data. The reason 
is that the asymmetries A{P^) are very sensitive to the xp integration range [x™™, x™'^^], and 
we cannot reproduce exactly the real [ ] of the data. This is supported by Fig. El 

in which we plot the calculated results of A{P^) with two options of XF-cut: the integrated 
ranges of [ [—0.12,0] and [0,0.12] respectively. We find that the magnitude 

of the asymmetries A{P^) differ significantly for the two options, whereas the trend of the 
P^-dependence is similar. Thus both the quark-diquark model and the pQCD based analysis 
can explain the P^-dependence in trend, but the exact magnitude is sensitive to the x^-cut. 
Therefore we conclude that the magnitude of A{P^) with detailed xp-cnt information is 
important for the purpose of distinguishing between different model predictions. Thus far 
we only expect our formalism to work for asymmetries, rather than for cross sections. In 
this last case there are questions of normalization and evolution which are not present for 
asymmetries. 

In summary, we showed in this work that the A-A asymmetries in hadron-nucleon coUi- 
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sions can provide information about the flavor structure of quark to A fragmentation func- 
tions. Both the quark-diquark model and a pQCD based analysis can explain the available 
data qualitatively. We also showed that the two models give significant different predictions 
for the magnitude of the Xi?-dependent asymmetries at large Xp for beams. Thus high 
precision measurements of the asymmetries with detailed Xp-cut and P^-cut information are 
important to discriminate between different predictions. 
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FIG. 1: The prediction for the A production asymmetries vs xp for various incident beams, and 
with three options for the contribution from unfavored fragmentation: (1) the dotted curve is 
for Cy = 1 and Cs = 0, a = 0;(2) the solid curve is for Cy = Cs = 1, a = 0.5; and (3) the 
dashed curve is for Cy = 1 and Cs = 3, a = 1. The left row figures correspond to results from 
the quark-diquark model, and the right row figures correspond to results from a pQCD based 
analysis. The integration range of is < < 3 (GeV/c)^ for positively charged beams, and 
< Pj' < 10 (GeV/c)^ for negatively charged beams. The experimental data are taken from 
Ref. 
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FIG. 2: The same as in Fig. ^ using only the contribution from unfavored fragmentation with 
Cv = C5 = 1 for a = 0.5. The integration range of is [0, 5] GeV^ for the thick sohd curves and 
[5, 10] GeV^ for the thin sohd curves. 



12 





(a) Diq. 

n* Beam 











(b) QCD 

n* Beam 










(d) QCD 

7t" Beam 





_1_ 




■ (g) 


Diq, 

K" Beom 











(i) Diq. 

p Beam 







(f) QCD 

K"^ Beam 















(h) 


QCD 

K' Beom 






FIG. 3: The prediction of A production asymmetries vs for various incident beams with the 
dotted, sohd, and dashed curves corresponding to the three options for the unfavored fragmentation 
as in Fig. n The integration range of xp is —0.12 < xp < 0.12 for positively charged beams, and 
—0.12 < Xp < 0.4 for negatively charged beams. The experimental data are taken from Ref. |^. 
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FIG. 4: The same as in Fig. IHl by using only the contribution from unfavored fragmentation with 
Cy = C5 = 1 for a = 0.5. The integration range oi xp is [-0.12, 0.0] for the thick sohd curves and 
[0.0, 0.12] for the thin soUd curves. 
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